arXiv:1506.03712vl [astro-ph.SR] llJun2015 


High resolution imaging of NGC 2346 with GSAOI/GeMS: disentangling the 
planetary nebula molecular structure to understand its origin and evolution. 

Arturo Manchado^’^’^, Letizia Stanghellini^, Eva Villaver^, Guillermo Garcia-Segura®, Richard A. 

Shaw‘S Sz D. A. Garcia-Hernandez^’^ 

ABSTRACT 

We present high spatial resolution (~ 60-90 milliarcseconds) images of the molecular 
hydrogen emission in the Planetary Nebula (PN) NGG 2346. The data were acquired 
during the System Verification of the Gemini Multi-Conjugate Adaptive Optics System 
-|- Gemini South Adaptive Optics Imager. At the distance of NGC 2346, 700 pc, the 
physical resolution corresponds to ~ 56 AU, which is slightly higher than that an [N II] 
image of NGC 2346 obtained with HST/WFPC2. With this unprecedented resolution 
we were able to study in detail the structure of the H 2 gas within the nebula for the 
first time. We found it to be composed of knots and filaments, which at lower resolution 
had appeared to be a uniform torus of material. We explain how the formation of the 
clumps and filaments in this PN is consistent with a mechanism in which a central 
hot bubble of nebular gas surrounding the central star has been depressurized, and the 
thermal pressure of the photoionized region drives the fragmentation of the swept-up 
shell. 

Subject headings: hydrodynamics-ISM: structure-ISM: jets and outflows-planetary neb¬ 
ulae: general-planetary nebulae: individual (NGC2346)-stars: AGB and post-AGB- 
stars: winds, outflows 


1. Introduction 


Molecular hydrogeii w as first detected in the Galactic planetary nebula (PN) NGC 
by IZuckerman fc Gatlevi (Il988l f . Thanks in part to its proximity, since then several authors 
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have 
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studi e d its Hp gas emiss ion in detail fe.g.. iKastner et alJll994l . Il996l : iLatter et alJll995l : IVicini et al 


I 999 I : lArias et alJ l200lh . They described an apparent structure of clumps and a torus not un¬ 


common in PNe. Neutral clumps have been detec ted in PNe in nearby extended objects such as 


NGC 6543 or NGC 6852 (jMeaburn &: LoDeall993ll and cometary knots seem to be det ected often 


in evolved PNe if the resolution allows it (e.g., iQ’Dell et alJl2002l : IHuggins et al.ll2002h . However, 


the formation of these small scale structures is not well understood and different mechanisms have 
been pr oposed to ex plain their formation. Some involve the formation of clumps prior to photoion¬ 


ization (jSokei:ill998l '). while others have suggested that pre-existing high density structures in the 


PN shell ( 

Aluzas et al. 

20121. The hvdrodvnamical evolution of those clumps has been studied bv 

Redman et al. ( 

20031 and 

Aluzas et al. 

(20121. and other mechanisms such as their formation in 


the p r ogenitor asymptotic g i ant b ranch (AGB) atmosphere has been ruled out (jHuggins &: Mauron 
2002l l. iGarcia-Segura et al.l (j2006l l proposed that clumps are the result of the fragmentation of the 
swept-up shell. 


1981 


hereafter 


The central star of NGG 2346 is in a close binary system (jMendez fc Niemelal 
MN81) with a measured period of ~ 16 days. Given its small orbital separation, it is possible that 
the system is the remnant of common-env elope (CE) evolution. Howeve r, alt ernative scenarios 
for its formation have also been proposed. iDavis. Kolb fc WillemsI ([20101), and Ide Kool fc Ritter 
(I 1993 I I suggested that the system did not evolve through CE, rather, it went through mass transfer 
from an evolved primary with a radiative envelope into the companion, via a thermally unstable 
Roche Lobe overflow. 


The role of binaries in the evolution from the AGB to the PN phase is not well under- 
stood. This is mainly because only 40 binary central-stars of Planetary Nebulae have been studied 
(jPe Marco et al.ll2013r). From these only a small fraction are beli eved to have undergone CE evolu¬ 


tion fe.g.. ICorradi et al.ll2014l . iMiszalski. Boffin &: Corradill2013lL In the CE phase the ejection o f 


the envelope is conhn ed closely to the equator (e.g., ISandauist et al.lll998l . iRicker fc Taamll2012lL 


Ricker fc TaamI (12012l ) found that 90% of the outflowing material is confined to an angle of 30° on 


either side of the equatorial plane. 


It this paper we study in detail the structure and gas distribution of H 2 in NGC 2346 using 
exquisite, high resolution images from observations with the Gemini Multi-Conjugate Adaptive 
Optics System (GeMS) -|- Gemini South Adaptive Optics Imager (GSAOI). We also study the 
properties of the binary central star system. In section §2 we describe the observations; Section §3 
includes a description of the molecular hydrogen emission; in §4 we discuss the stellar progenitor 
mass; section §5 gives a general discussion of the results; hnally, we summarize our hndings and 
conclusions in S6. 
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2. Observations 


We used the GeMS + GSAOI during System Verification (March 2013). GeMS (jd’Orgeville et ah 


2012 l i comprises multiple deformable mirrors, and uses three Natural Guide Stars (NGS) and five 


sodium Laser Guide Stars (LGS) . GSAOI is a near infrared (NIR) camera used with GeMS on 
Gemini South (jPessev et ahl 120131 ). GSAOI provides diffraction limited images in the 0.9—2.4 /rm 
range, using a 2 x 2 mosaic Rockwell HAWAII-2RG 2048 x 2048 arrays. The GSAOI field of view 
(FOV) is 85'' X 85" with a scale of 0.02" pixel“^ and a gap between the arrays of ~ 2mm (which 
corresponds to 2 . 2 "on the sky). 


We obtained narrow-band images in the H 2 (1-0) S(l) 2.122 /rm, Br 7 2.166 /rm, and H 2 (2-1) S(l) 
2.248 /rm filters with good natural seeing conditions (0.4-0.5"). We used 5 dither positions with indi¬ 
vidual exposure times of 120 s for the H 2 (1-0) S(l) images and 360 s for the Br 7 and H 2 (2-1) S(l) 
images. As the object is larger than the FOV, adjacent sky frames were taken with the same ex¬ 
posure time. Exposure time, on-target, was 600 s for the H 2 (1-0) S(l) image, and 1800 s for the 
Br 7 and H 2 (2-1) S(l) images. Data reduction, including distortion correction, was carried out 
using the Gemini IRAF0 package vl.l2, gsaoi. After the combination of all the images for each 
band we found the FWHM varies over the whole FOV, between 60 and 90 mas, with an average 
value of 80 mas. We used 2MASS band magnitudes of several field stars to flux calibrate the 
H 2 (TO) S(l) images. 


Transmission of the H 2 (1-0) S(l) 2.122 /rm filter is about 90% and the width is 0.032 /rm. A 
color correction between the center of the band (2.159 /rm), and the center of the H 2 (1-0) S(l) 
line (2.122 ^m), line was applied, which led to an uncertainty of 1% in the flux calibration. The 
statistical error in the flux calibration, as measured from different field stars is about 2%. Thus, al¬ 
together, the error in the flux calibration is approximately 2.2%. We obtained an inverse sensitivity 
for point-sources of 4.94 ±0.10 x 10“^^ erg“^ cm“^ s“^ per electron. The pixel area is 9.5 x 10“^® sr, 
so the inverse sensitivity for extended sources is 5.25 ±0.10 x 10““^ erg“^ cm“^ s“^ sr“^ per electron. 


Fig. [T] shows the resulting image of the H 2 (1-0) S(l) 2.122 /rm emission. The signal-to-noise 
(S/N) ratio of the extended emission in the Br 7 and the H 2 (2-1) S(l) images is less than 2a and 
they are not shown along the paper. The [N II] image was obtained with the WFPG2 on board the 
HST (Proposal SM^ERO-7129). The images were retrieved from the Mikulski Archive for Space 
Telescopes (MAST)g. In dividual images were combined in order to remove the cosmic rays using 
the DrizzlePac software ([Gonzaga et al.ll2012l L with a resulting exposure time of 1280 s. 


^linage Reduction and Analysis Facility (IRAF) software is distributed by the National Optical Astronomy Ob¬ 
servatories, which is operated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 

^https://ar chive.stsci.edu/hst/ 
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3. H 2 emission 


was first detected in 

NGC 2346 bv 

Zuckerman &: Gatlev 

(19881 

Kastner et al. 

(1994 


Vicini et al. 

(. 

1999 

) and 

Arias et al. 

(20011 


using a circular variable filter, 
imaged the H 2 (1-0) S(l) emission in the whole nebula, with a seeing of about 2"5. They resolved 
the extended H 2 (1-0) S(l) emission in both the lobes and what they thought was a central torus. 
However our much higher resolution ima ges show tha t this feature is not a smooth torus, but an 
aggregation of a large number of clumps. IVicini et al.l (jl999l ) also obtained spectra in the K-band, 
obtaining a rati o of the H 9 (TO) S fD/H? ( 2 - 1 ) S(l) lines of about 14 in the central part and 4.3 in 
the lobes, while lArias et al.l (| 200 ll ) found a value of 6.6 for the torus. Our much better resolution 
images (see Fig. show that the H 2 (1-0) S(l) emission is in the form of clumps and cometary 
knots (see Fig. [2] and Fig. [3|), and are located in both the lobes and the central region. Clump sizes 
vary from Cf.'16 to 0l'34 along the semi-major axis, and projected distances from the central star are 
from 1"6 to 44"2. For the brightest clump we derive a H 2 (1-0) S (l)/H 2 (2-1) S(l) line ratio of 4.7, 


which is consistent with the value derived bv IVicini et al.l (jl999l l for the lobes and bv lArias et al 


(j200ll ) for the central region. From our H 2 (1-0) S(l) image we obtain line surface brightness of 


8.8 X 10 ^ and 1.8 x 10 ^ erg ^ cm ^ s ^ sr ^ for the brig htest and weakest clumps, respectively. 


-1-4 


v-i, 


,-2 




The highest value is 9 times higher than IVicini et al.l (j 19991 ) found. This is due to the dilutio n effect 
of their poor spatial resolution. If we degrade our H 2 (1-0) S(l) image to lVicini et all (j 19991 ) spatial 
resolution (3^^) we obtain a very similar value of 0.9 x 10“^ erg“^ cm“^ s“^ sr“^ compared with 
their value of 0.7 ± 2 x 10““^ erg“^ cm“^s“^sr“^. The val ues that we hnd fo r the brightest clumps 
are 2.5 times higher than the predictions of Mod el 1 fromIVicini et al.l (jl999l l in the 3000-10000 yr 
in terval, and 5 times higher than the value that ISoeck et al.l (j2003l ) found in NGC 6720. Model 1 


m 


Vicini et al.l (jl999l l predicts H 2 (1-0) S(l) emission originating in the photodissociation region 


(PDR), created at the edge of the neutra l shell by the UV radiation. Thus, our values seems to favor 
shock excitation. However IVicini et al.l (jl999l ) model reproduced the average surface brightness, 
and did not explore the full parameter space. 


4. Stellar Progenitor Mass 


The central star of NGC 2346 is a spectroscopic binary (MN81) with a period of 15.99 days. 
The main sequence companion of the ionizing star, which is seen in Fig. (TJ is a 1.8 Mq A5V star 
with Tg// = 8000 K and L = 18 Lq (MN81). From the magnitude and extinction of the companion, 
MN81 obtained a distance of 700 pc to the system. The ultraviolet excess of the ionizing star is 
compatible with a white dwarf with L = 50 Lq and Tg/j = 100,000 K. 


The mass of the ionizing star carries the uncertainty in the unknown inclination angle of the 
binary orbit, but a reasonable assumption is to consider the orbital plane to be perpendicular to 


the bipolar lobes of the nebula and zero orbital eccentricity (| Jones et al 


inclination angle of the lobes with respect to the line of sight estimated by 


20inlL Given the 120° 


Arias et ahl ( 200 l|), the 
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binary orbit al plane would be seen at an angle of i = 60° and within 35° < i < 85° using a 25° 
uncertainty (j Arias et alJ 1200 li b In Tabled] we list the solutions for the mass of the ionizing source 
and the separation of the two components using this plausible range of inclination angles for the 
binary system. 

Using the plausible range of inclination angles, we obtain a mass range for the ionizing source 
between 0.32 and 0.72 Mq. The binary separation, which has a very small dependency on the 
inclination assumed, has values ranging from 0.159 to 0.167 AU (34.2 — 35.7 R©). The ionizing 
star must have evolved past the giant phase, reaching a radius (either on the red giant or on the 
AGB) much larger than the current orbital separation of the system. Thus, it is very likely that 
the binary system experienced orbital decay if it suffered CE evolution, but significant orbital 
decay is not expected if instead mass was simply transferred from the primary when the giant’s 
envelope was not fully convective. The evolution of a binary system through the CE phase is not 
und erstood in detail, an d thus it is difficult to infer the precise history of individual systems (see 
e.g., Ilvanova et al.ll2ni3l ). Given the orbital parameters of the known components and the mass of 


the main sequence companion to the PN progenitor, there are two reasonable assumptions we can 
make regarding the past history of this system: i) the evolution of the binary did not result in a 
merger; and ii) the main sequence mass of the PN progenitor must have been greater than 1.8 M© 
for it to have evolved into a giant before its companion. 

It is difficult to estimate the amount of mass accumulated by the secondary. Note that the 
fraction of the envelope that was lost in the CE evolution is unknown. Therefore we can only 
set a lower limit to the initial mass of the PN progenitor. Likewise, we have to assume that the 
initial orbit separation was larger than 0.16 AU. But since the spectrum of the ionizing source is 
unknown, we are unable to determine with certainty whether the CE phase occurred during the 
progenitor’s red giant branch (RGB) or AGB evolution. However, we argue that the CE phase 
must have occurred on th e AGB, since the orbital separation is expected to decrease during CE 
evolution (jWebbinkI 1 19841 1 an d considering that the largest core mass that can build up during the 
RGB is ~ 0.47 M© (see e.g. ISweigart et al.l (Il990l )). This would give the progenitor enough time 
to build up a large enough CO core to sustain the photoionization of the nebula. 


Davis. Kolb &: WillemsI (j2010ll . using population synthesis techniques, modeled a CE system 


to obtain ZAMS progenitor masses of 2.47 M© and 0.98 M© and a final WD configuration with 
mass 0.33 M©. Such a small WD mass however cannot sustain the ionization of the nebula, and it is 
not consistent with the UV emission observed. A more m assive primary tha t expe rienced enhanced 
mass-loss before it filled its Roche lobe was suggested by iTout fc EggletonI (jl988l l for this system. 


Since the companion has mass of 1.8 M© , we c an set this value as the lower limit to the 
PN progenitor mass. Erom IVassiliadis &: Wood! (|1993l l this value is consistent with L = 50 © and 
Tg// = 100,000 K. Thus, the most plausible solution for the progenitor is a 1.8 M© star that has 
experienced CE evolution on the AGB, leaving a remnant greater than 0.61 M©. 


We computed the number of ionizing photons for a 1.8 M© star assuming a black-body spec- 























trum, log L/Lq = 1.730, and Te// = 93,540 K, to be 4.0 x 10^^ photons s We will use this 
result in the modeling described in the next section. 


5. Discussion 


5.1. Analysis of the Clumps 


As can be seen in Figures [T] through [H the H 2 emission is clumpy. This clumpy structure is 
in contrast, and as expected, with the smoother structure seen in the ionized gas (Fig. [1]). Our 
high resolution images clearly show that the H 2 emission is not uniformly distributed in a torus 
around the central star, as it has been interpreted from previous, lower resolution observations, 
but displays a fragmented structure mainly composed of clumps and cometary knots. To further 
probe this point we convolved our H 2 (1-0) S(l) image with a symmetric PSF of 2", as shown in 
Fig. [5l It is easy to see how the presence of a torus around the central star might be inferred from 
a lower resolution image. Most of the clumps are generally superimposed on a diffuse emission (see 
Figures [T] through H]). They are sometimes grouped as in Fig. [21 have cometary tails as in Fig. [3l 
or have bow shock shape as in Fig. |4l All of this makes it very difficult to measure the number and 
distribution of clumps in the nebula (e.g., total number of clumps, radial distribution, fluxes and 
knot structure). In any cas e, the H 2 (1-0) S(l) d istribution in NGC 2346 is completely different 
than the one in NGC 7293 ( Matsuura et al. 2009 ) 


Assuming a distance of 700 pc, each pixel in our image corresponds to a physical size of 14 AU. 
From Fig. [T]we find that most of the knots are concentrated on the waist of the nebula (East-West), 
and at the edges of the lobes. In the lobes, at distances further than 0.02 pc from the central star, 
the H 2 emission is more diffuse. The clumps are at distances of 0.0089 to 0.15 pc from the central 
star. Typical sizes are in the 1.67-3.52 x 10^® cm (112 to 238 AU) range. 


5.2. Formation of the Clumps 


Several mechanisms have been explored to explain the formation of molecular clumps and 
cometary tails in PNe. A scenario in w hich pre-existing h igh density structures in the ISM are 
responsible for the structures observed (jAluzas et al.ll2012h is difficult to rule out a priori. The 
question is t hen how the ISM high den sity clumps were formed in the first place. The mechanism 
proposed bv lGarcia-Segura et al.l (j2006l l provides a solution for both the formation of the observed 
H 2 clumps in the PN and the possible existence of long-life high density clumps in the ISM. In 
addition, w hen high resolution observ ations of H 2 in PNe are obtained, they often reveal clumpy 
structures (jMarquez-Lugo et al.ll2r)13l ). To assume that the ejected AGB shell would always en¬ 
counter a clumpy ISM does not seem reasonable; it would imply the ISM to be densely populated by 
high density structures along many lines of sight. We prefer the more plausible explanation offered 
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by iGarcfa-Seeura et alJ (jl999l ) that Rayleigh-Taylor instabilities, or the effects o f the ionizing radi¬ 
ation in the nebula, lead to the formation of clumps. iGarci'a-Seeura et al.l (j2006l i proposed that the 
cessation of a fast stellar wind can self-consistently produce the type of fragmented structures we 
observed in NGG 2346. The mechanism only requires a rapid decline (i.e., switching off) of the fast 
stellar wind in an ionization-bounded PN. If the stellar wind becomes negligible the hot, shocked 
bubble depressurizes, and the thermal pressure of the photoionized region, at the inner edge of 
the swept-up shell, becomes dominant. The shell tends to fragment, creating clumps of neutral 
gas with comet-like tails and long, photoionized trails in between, while the p hotoionized material 


expan ds back toward the central star. In the hydrodynamical simulations by iGarcfa-Seeura et al 


(|2006l i the cometary globules originate in the neutral, swept-up shell of piled up AGB wind and so 
they must be excited by shocks. 

There are two mechanisms arising in very diffe rent physical conditions th at can be responsible 
for the observed H 2 emission: (i) s hock excitation (IShull &: Hollenbachlll978l) . or (h ) UV pumping 
in photo-dissociation regions (e.g.. iBlack fc van Dishoecld 119871 ). I Arias et al.l (120011 ) have explored 
the nature of the excitation in NGG 2346 using the relative intensities of the H 2 (1-0) S(l) and 
the H 2 (2-1) S(l) lines, and concluded that the H 2 emission line ratios are consistent with shock 
excitation, both in the central nebular regions and in the lobes. Further support for shock excitation 
comes from the incompatibility of H 2 (1-0) S(l) flux with UV excitation, described in §3. 

Since the spectrum of the central star is unknown, we cannot determine for certain if the 
antecedent fast wind has decayed, leading to depressurization of the hot bubble. However, there 
is some evidence of this being the c ase. First, the cent ral star is relatively massive, so the wind 
would be expected to evolve quickly (jVillaver et al.ll2002lb Second, in spite of moderate interstellar 
extinction, if a robust, fast wind were present on e would expect X-ray emission. Yet NGG 2346 
has not been detecte d with either XMM imaging (jGruendl et al.ll2006l i nor Chandra observations 
(jKastner et al.ll2012l b Furthermore, orbita l decay is compatible with the nebula kin ematical age 
(> 3,500 yr) derived bv I Arias et al.l ([20011), which, according to IVillaver et al.l (j2002lb is evidence 
that the fast wind has already declined below its maximum kinetic energy. 


T o further examine the formation of the clumps, we have used the 3D version of the lGarcfa-Segura et al 
(2006) model, with AGB wind velocity Vexp = 10 kms“^, AGB mass loss M = 10“® Mq yr“^, 
fast wind velocity Vexp = 1000 kms“^, and fast wind mass loss M = 10“'^ Mq yr“^. The tem¬ 
poral grid has 125^ zones, and physical dimensions 0.8 x 0.8 x 0.8 pc. The fast wind was termi- 
nated in the simulation at 1000 yr from the onset. The photoionization follows the prescription of 
Garcfa-Segura fc Francol (199a)) assuming a central ionizing emission of lO'^® s“^ ionizing photons 
(see above). To obtain bipolarity, we assurned tha t the rotation velocity approached 95% of its 
critical value, following iGarcfa-Segura et al.l (119991 ). The simulation was done using the hydro- 


dynamical code ZEUS-3D (iNormanI 120011 ) . It must be noted here that this simulation is highly 


qualitative, due to the following factors: one, the grid resolution was low, and two, the simulation 
was not tailored to explain NGG 2346 in detail but instead used standard average parameters (al¬ 
though the model was optimized for NGG 2346 in that it assumed a similar number of phoionizing 























































photons as those derived from the NGC 2346’s UV flux). Nevertheless, this model is useful for 
understanding clump formation. 

In Fig. [6] we show the gas emission measure for nine model snapshots, on the X-Z plane. The 
initial model is for t=1000 yr (following the termination of the fast stellar wind), and subsequent 
models are at timestep of 1,000 yr, covering 9,000 yr of the evolution. The emission of the various 
components are represented by different colors: photoionized gas with temperature~10,000 K is 
represented in green; gas (either photo ionized or neutral) with temperature above 1,000 K is in red, 
and gas with temperature below 1,000 K (blue). The blue (cool) gas emission follows the molecular 
emission (i.e., H 2 ), so we can use this latter indicator to see what happens in the molecular regime. 
The simulation shows, as a general view, an initial toroid of cool gas at the equator. Once the 
swept-up shell is highly fragmented, the toroid is no longer visible and only the large clumps, with 
an optical depth big enough to shield the ionizing radiation, will survive and be detected. 


The current evolutionary stage of NGC 2346 is represented by models of times between the 
hrst and the second snapshot of Fig. [6l since it appears that the clumps in NGC 2346 are just being 
formed, and the photoionization front has not yet reached the equatorial latitudes. NGC 2346 has 
a small number of long cometary knots (Fig. [3]), not easily recognable in the HST [N ii] image 
(Fig. d]). This contrasts with the large number of such knots in the Helix (NGC 7293) or in the 
Eskimo (NGC 2392) nebulae ( O’Dell et al. 2002 1. In addition the siz e of the clumps in NGC 7293 
(94 AU) are similar in size to those in NGC 2346 ( O’Dell et al. 2007 ). 


5.3. The Survival of H 2 Clumps 


An interesting question is to determine how the molecular hydrogen clumps can survive the 
photoionization front. If the propagation velocity dR/dt of the ioniz ation front within a clump can 
be expressed according to equation 7.5 of iDvson fc WilliamsI (1199711 . then 


dR S, 1 ^ ^ 

TTT = - o -^^0P2, 

dt 47r no 3 


( 1 ) 


where S'* is the number of photoionizing photons, R the distance from the ionizing source to the 
clump, dR the size of the clump, /32 the recombination coefficient, and no the density of the clumps. 
For NGC 2346, the closest clumps (see Fig. |3|) are at a projected distance of 2.75 x 10^® cm (1850 
AU) from the central star, with size ~ 2.1 x 10^^ cm. Thus, assuming S* = 4.0 x 10^® s“^ and the 
kinematical age as a proxy for the evolutionary time (3500 yr), we estimate the minimum density 
required for clumps to survive the passage of the ionization front. We obtain a value of no = 15, 000 
cm“^, which is smaller than that of the gas density derived for the largest clumps (5 x 10^ cm 


- 3 ^ 


from dust absorption of the [O III] images in NGC 6853 bv iMeabnrn fc Lopea (119931 ). 


During post-AGB evolution the stellar luminosity declines, and some nebular clumps may not 
be photoionized. Instead, they survive as molecular material, eventually making the ISM clumpy. 
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Since the H 2 (1-0) S(l) will not couple with the 3 K cosmic microwave background, the only 


destr uc tion mechanism in the ISM is the spon taneous radiative dissociation (jStephens &: Dalearno 


197l|). iPfenniger. Combes &: Martinetl (|l994l l proposed that H 2 clumps of a size of about 30 AU 
can explain baryonic dark matter in spiral galaxies. The H 2 clumps in NGC 2346, with sizes 
are between 112-238 AU, will populate the ISM, contributing to the baryonic dark matter of the 
Galaxy. Whether H 2 (1-0) S(l) in clumps survive or not past the PN stage and populate the ISM 
of galaxies is beyond the scope of this work. We speculate that if clumps survive, given the amount 
of mass they contain, they should be accounted for as a source of baryonic dark matter. 


6. Conclusions 

We have obtained molecular hydrogen (H 2 ) images of NGG 2346 with unprecedented spatial 
resolution; 56 AU at an adopted distance of 700 pc. The images reveal the H 2 emission to be 
fragmented in clumps and cometary knots, rather than a uniform disk as previously thought. The 
clumps range in size from 112 to 238 AU; the clump apparently closest to the central star lies at a 
projected distance of 1850 AU. The central star has undergone binary interaction with its nearby 
companion, probably during its AGB phase. The ionizing star has a minimum mass of 0.61 Mq, 
and an ionizing flux of 4.0 x 10^® photons s“^. 

We performed a hydrodynamical simulation that shows how an initial disk or toroid breaks 
up into individual clumps, once the swept-up shell fragments following the decline of the fast 
wind. In order to survive the ionization front, the pre-ionization density of the molecular hydrogen 
clumps must exceed ~ 15,000 cm“^, which is found in some PNe. Those clumps that survive the 
ionization front will eventually populate the ISM, and may contribute to the baryonic dark matter 
in our Galaxy. 
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Table 1. Masses and separation of the binary system 


Inclination Angle 
(deg) 

Progenitor Mass 
(Mo) 

Separation 

(AU) 

30 

0.72 

0.167 

35 

0.61 

0.166 

38 

0.56 

0.165 

40 

0.53 

0.164 

50 

0.43 

0.162 

60 

0.37 

0.161 

70 

0.34 

0.160 

80 

0.33 

0.159 

85 

0.32 

0.159 
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Fig. 1.— Images of NGC 2346 obtained with GSAOI in H 2 (1-0) S(l) {left panel), and with 
HST/WFPC2 in [N ii] (right). Images are oriented with N up and E to the left; tick labels are in 
arcseconds with respect to the central star. 
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Fig. 2.— H 2 (1-0) S(l) clumps NW of the central star. 
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Fig. 3.— H 2 (1-0) S(l) cometary knot NW of the central star. 
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Fig. 6.— Snapshots of the emission measure (cm“® pc) along the line sight of the modelled gas. 
Three different temperatures are represented in different colors (see text for details). Starting at 
1000 yr elapsed time with the cessation of the stellar wind {upper left), the evolution proceeds from 
upper to lower, then left-to-right, with a timestep of 1000 yr. 



